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Pade Unitarizations: A Critical Look 
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We use the Linear Sigma Model to test the convergence of different sequences of Pade Approximants when used 
to unitarize the low energy mr scattering amplitude. We find that, in this particular case, diagonal sequences 
reproduce the sigma pole with high accuracy, as opposed to other sequences that have been discussed extensively. 
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1. INTRODUCTION 

One of the most interesting aspects of strong 
interactions is the enormously rich and compli- 
cated hadronic spectrum that arises from such 
a simple high energy structure. How does ex- 
actly this spectrum arise from QCD is a ques- 
tion under study. The understanding of QCD at 
low energies is basically provided by ChPT [I], 
the low energy effective theory of QCD, which 
gives also some insight on the spectrum issue. 
Although ChPT breaks down in the vicinity of 
the resonance region, arguments following from 
unitarity and analyticity are able to extract in- 
formation on some resonances from the chiral La- 
grangian, which contains such information hidden 
in the values of its low energy constants. One of 
the most successful of such unitarization proce- 
dures is the Inverse Amplitude Method (IAM) [5] , 
which is able to extract with acceptable accuracy 
the lowest lying scalar and vector resonances. 

These proceedings are based on the work of 
Ref. 3], where some of the unitarization tech- 
niques are tested by means of the Linear Sigma 
Model, for which exact information about the 
sigma resonance is available. We focus on the 
IAM and also on some sequences of Pade Approx- 
imants (PAs) [4]. 

2. UNITARIZATION OF THE NLSM 

In order to determine the scalar meson mass 
and width up to O(g), we compute the one-loop 
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where p(s) = y/l — AM 2 / s , and the term —13/3 
is determined by the renormalization scheme that 
sets the relation 2gF 2 — M 2 at the one-loop or- 
der, with F the pion decay constant. Now it is 
possible to extract the pole s p of the propagator 
up to the considered order in perturbation the- 
ory. If one approaches the branch cut from the 
upper part of the complex s-plane, the pole in 
the second Riemann sheet is located at 
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up to corrections of (D(g 2 ). The pole mass and 
width, defined from s p — (M p — iT p /2) 2 , are then 
given by 
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We now consider the LSM at low energies. The 
contribution from the sigma exchanges to the 
rcnormalized 0(p ) couplings gives [1] 
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The 7T7r-scattering is determined by the tt + tt~ — > 
7r°7r n amplitude, which is given up to 0(p 4 ) in the 
chiral expansion by 
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where /i refers here to the arbitrary renormaliza- 
tion scale, and the chiral limit has been consid- 
ered. With this one constructs the definite isospin 
partial waves, tj(s), with I J = 00, 11, 20, finding 
the following 0(p 2 ) amplitudes, 
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up to corrections of 0(g 2 ). Given the 0(p 2 ) 
and 0(p 4 ) xPT amplitudes from Eqs. 
it is then possible to extract the poles of the 
corresponding t(s) IAM for the LSM, satisfying 
i ( s )(4)A(s)(2) = 1 at s 
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These are the poles that appear in the unphysical 
Ricmann sheet as one approaches from upper half 
of the first Riemann sheet. There is also a conju- 
gate pole at s* if one approaches the real s-axis 
from below. 

The first thing to be noticed is that poles ap- 
pear in the I J = 11 and 20 channels even for 
small values of g, contrary to what one expects 
in the LSM, where no meson with these quan- 
tum numbers exists. Furthermore, these "states" 
are not resonances, as they are located on the 
left-hand side of the complex s-plane, out of the 
physical Riemann sheet, and carrying a negative 
squared mass. 

As for the I, J = 00 channel, one finds a reso- 
nance with pole mass and width, 
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The I AM predictions for M 2 and M p T p result, re- 
spectively, 40% and 80% smaller than the original 
ones m the LSM, computed in Eq. (JSJ). There- 
fore, this provides an example of a situation in 
which this particular method cannot reproduce 
the hadronic properties of the theory from its ef- 
fective low-energy description. 

3. HIGHER ORDER PA's 

Now we consider higher order PAs. In order to 
be able to handle the amplitude at higher orders, 
we will consider the irir scattering at tree-level. 
This is equivalent to working in the limit j« 1 
and keeping just the first non-trivial contribution 
in the g expansion. At low energies the amplitude 
becomes 
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Figure 1. Position of the nearest pole to M 2 for 
the first PAs of the form [l/N] with A odd (for 
even A all the poles are complex). 

and the partial waves are given by, 
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In this section we compare the [l/N] and [N/N] 
sequences for the study of the tttt partial wave 
scattering amplitudes. We shall focus on the 
IJ — 00 partial wave, but analogous results 
are found for the other channels. Former works 
pointed out that the PAs and other unitarizations 
fail to incorporate the crossed channel resonance 
exchanges. Nonetheless, we will see that as N 
grows, the poles of the sequence [N/N] actually 
tend to mimic not only the s-channel poles but 
also the left-hand cut contribution from diagrams 
with resonances in the t and u channels. 

We begin with the sequences of the type [1, A]. 
Our results are summarized in Figs. [1] and [5] No 
convergence is found with this sequence. In the 
case of N odd, Fig. Q] shows that the P^ pole 
closest to M 2 does not approach this value even 
for very large A, always remaining a 30% below. 
The analytical structure of the original amplitude 
(s-channel sigma pole plus left-hand cut) is never 
recovered since the [1 /A] PAs always set the poles 
in the circular pattern shown in Fig. O This sug- 
gests that the use of further [1/A] approximants 



Figure 2. Poles of the P^ in the complex plane. 

to extend the IAM is not the optimal way to pro- 
ceed, even if we had an accurate knowledge of the 
low-energy expansion up to very high orders. 

Alternatively, the use of sequences such as 
[N + K/N] (e.g. [A-2/A], [A-l/A], [A/A]...) 
seems to be a better strategy. In the following we 
analyze the sequence [N/N] , as it ensures the ap- 
propriate behavior at high energies, \t(s)\ < 1. 
Nevertheless, similar results have been generally 
found for the [A + A/A] PAs with K ^ 0. The 
P$ pole closest to M 2 is shown in Fig. 03 One 
finds a quick convergence of the sequence: Pi 
reproduces the sigma pole a 40% off but P 2 dis- 
agrees by less than 1%, P| by less than 0.1%, 
etc. Notice that already P 2 2 provides a much bet- 
ter description than Pg 1; although one includes 
far more low-energy information in the latter. 
All this points out the sizable discrepancy of the 
first element of the sequence (Pi) with respect to 
the original amplitude. It also indicates that the 
[1/A] PAs do not produce a serious improvement. 
On the contrary, the [A/A] sequences provide a 
more efficient strategy with a quick convergence. 

Likewise, Fig. [5] shows how the [1/A] PAs are 
unable to recover the analytical structure of the 
original amplitude, whereas the [N/N] sequence, 
besides providing the isolated pole of the sigma, 
tends to reproduce the left-hand cut as A in- 
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are plotted in Fig. 2J 
Although a PA is a rational function without cuts, 
these are mimicked by placing poles where the 
cuts should lie. The P 2 2 o nas one isolated pole 
near M 2 (with an accuracy of 10 -30 ) and nine- 
teen poles over the real axis at s p < —M 2 , i.e. 
on the left-hand cut of the original function. As 



Figure 3. Location of the closest pole to for 
the first [N/N] PAs. 

iV is increased, the number of poles lying on the 
branch cut increases too. 

As an amusement, we have also probed the PA 
sequence [N / 1] which has the same number of in- 
puts as the [1/iV] has for a given N. In this new 
case we have found convergence in both LSM and 
the resonance model presented in the following 
section but slower than the [N/N]. For instance, 
the prediction for the M% for the first are 
j^r = 0.55, 1.47, 0.73, 1.27, 0.81... A criticism that 
can be done to this sequence is its lack of unitar- 
ity, in contrast to the other studied sequences. 

4. CONCLUSIONS 

The conclusion seems to be two-fold. First, the 
LSM features a low energy behavior for which the 
IAM does not provide a good approximation, at 
least in what concerns the sigma pole. Second, 
while some PAs converge amazingly well, even be- 
yond expectation (there are in these cases no the- 
orems ensuring convergence), as for example the 
diagonal sequences [N, N] , some other sequences 
such as [1 , N] show a particularly poor conver- 
gence. In addition, the convergent sequences dis- 
play an analytic structure that mimics perfectly 
poles and cuts. 

Concerning the first point, some explanation 
has been provided by some of the authors of [2] . 
They argue that the IAM is a good approxima- 
tion in strongly interacting regimes where unitar- 
ity violations are relevant, while in the case of 
a weakly interacting set-up as the one described 
here, in the limit g -C 1 the chiral expansion 
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Figure 4. Poles of P^q. 

breaks down. This is due to the fact that in this 
limit the contact C(p 4 ) terms are as important as 
the loop contributions. The IAM would behave 
properly for values of the LECS of order 10 -3 , as 
is the case in QCD. Moreover, they argue that in 
the weakly interacting limit the sigma resonance 
is almost decoupled, and thus the fact that the 
IAM can pin it down with a precision of about 
40% is remarkable. A more careful investigation 
of how the approximations involved in the deriva- 
tion of the IAM, break down in the LSM, would 
help to clarify this issue. 
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